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Abstract 

Context. X-ray binaries are usually divided in persistent and transient sources. For ultracompact X-ray binaries (UCXBs), the mass 
transfer rate is expected to be a strong function of orbital period, predicting persistent sources at short periods and transients at long 
periods. 

Aims. For 14 UCXBs including two candidates, we investigate the long-term variability and average bolometric luminosity with 
the purpose of learning how often a source can be expected to be visible above a given luminosity, and we compare the derived 
luminosities with the theoretical predictions. 

Methods. We use data from the RXTE All-Sky Monitor because of its long-term, unbiased observations. Many UCXBs are faint, i.e., 
they have a count rate at the noise level for most of the time. Still, information can be extracted from the data, either by using only 
reliable data points or by combining the bright-end variability behavior with the time-averaged luminosity. 

Results. Luminosity probability distributions show the fraction of time that a source emits above a given luminosity. All UCXBs 
show significant variability and relatively similar behavior, though the time-averaged luminosity implies higher variability in systems 
with an orbital period longer than 40 min. 

Conclusions. There is no large difference in the statistical luminosity behavior of what we usually call persistent and transient 
sources. UCXBs with an orbital period below ~ 30 min have a time-averaged bolometric luminosity that is in reasonable agreement 
with estimates based on the theoretical mass transfer rate. Around 40 min the lower bound on the time-averaged luminosity is similar 
to the luminosity based on the theoretical mass transfer rate, suggesting these sources are indeed faint when not detected. Above 50 
min some systems are much brighter than the theoretical mass transfer rate predicts, unless these systems have helium burning donors 
or lose additional angular momentum. 
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1. Introduction 

Ultracompact X-ray binaries (UCXBs) are low-mass X-ray bi- 
naries characterized by an orbital period of less than one hour. 
Given that the donor stars in these binary systems fill their Roche 
lobe, they must be white dwarfs or helium burning stars, since 
only these types of stars have an average density that corre- 



sponds to such a short orbital period (|Paczynski & Sienkiewicz 



|1981j [Nelson et al.|1986). The accretors can be neutron stars or 



black holes, though the latter have not been identified yet. 

Historically, X-ray sources including UCXBs have been cat- 
egorized as persistent or transient, based on whether a source 
was permanently visible by a given instrument or only occa- 
sionally. However, such a subdivision is an oversimplification, as 
sources have been found to vary over large ranges of timescales 
and amplitudes ( Levrne et al.|201 1) . One type of variability, first 
observed in dwarf novae, is attributed to a thermal-viscous in- 
stability in the accretion disk ( |Osaki||1974| |Smak||1984| |Lasota| 
|2001| i. As the disk grows in mass and heats up to close to 
the ionization temperature of the dominant element, the opacity 
strongly increases and radiation is trapped, and the disk's surface 
density approaches a threshold value. The temperature continues 
to rise more rapidly and the disk enters the hot state. The re- 
sulting high viscosity enhances the outward angular momentum 
transport in the disk, allowing for a higher accretion rate. This 
can be observed as an outburst, after which the cycle repeats. For 
UCXBs the mass transfer rate (and thus the disk temperature) is 



expected to decrease as a function of orbital period, suggesting 
that systems with orbital periods above ~ 30 min should undergo 
this instability in their outer disks (see Deloye & Bildsten 2003 ). 



1.1. Present research 

We investigate the bolometric luminosity distributions of the 
known UCXB population including candidates with tentative or- 
bital periods, which show the fraction of time a source spends at 
a given luminosity. The relation between the luminosity distri- 
bution and the average luminosity may give insight in the way 
transient and persistent behavior should be interpreted. 

We do not investigate periodicities in the light curves since 
this has been done already ( |Levine et al. 201 l| l, and periodicities 
are not relevant for the present purposes. 



2. Method 

Here we describe the known UCXB population and their obser- 
vational records, and how we use these data. 



2.1. Known UCXB population 

UCXBs are very rare objects. Presently 13 have been identified 
with high confidence (by finding the orbital period), all Galactic. 
These are located in globular clusters, the Galactic Bulge and the 
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Galactic Plane. Table[T]presents an updated version of the UCXB 
and UCXB candidat e overviews in|in'tZand et aL] ( |2007| ), |Lasota| 
|etaL| ( |2008) and |Nelemans & Jonker| l |2010) . 



2.2. Observations 

We use observations by the Rossi X-ray Timing Explorer All-Sky 
Monitor (RXTE ASM) ( |Bradt et al.|[T993) |Levine et al.1|T996) . 
The ASM has collected 16 yr X-ray Ught curves and there- 
fore gives insight in the source behavior over a large range of 
timescales. By its very nature of being an all-sky monitor, the 
ASM is an unbiased sampler of light curves, i.e., its observa- 
tions are not driven by source behavior. Pointed observations by 
e.g. Chandra or RXTE PCA have a higher sensitivity than the 
ASM, but are available far less regularly and often during atypi- 
cal source behavior. Hence, these observations cannot be used to 
gain insight in the 'normal' behavior over timescales of years. 



2.2.1. RXTE All-Sky Monitor 

The RXTE ASM has monitored the sky in the 2 - 10 keV en- 
ergy range from January 5, 1996 to January 3, 2012 (MJD 
50087 - 55929). Individual exposures, called dwells, have an 
integration time of typically 90 s. The intervals between dwells 
are irregular and on average ~ 2 hr. Light curves have been pub- 
lished in the ASM Products DatabasaMfor many resolved X-ray 
sources, including 12 of the 13 confirmed UCXBs. X-ray pho- 
ton count rates in the ASM energy range have been converted 
to bolometric flux using the c onversion of 1 ASM cou nts s = 
7.7 x 10" 10 erg cirT 2 s" 1 by |in't Zand et al.| ( |2007| ) (with an 
error of a factor of ~ 2) and to bolometric luminosities using 
the distance estimates in Table [Tj The ASM does n ot provide a 
light curve for NGC 6440 X-2 ( |Heinke et al.pjlO) and UCXB 
candidate NGC 6652 B. Both share their host globular cluster 
with a brighter X-ray source. The latter, NGC 6652 B, is too 
close to XB 1832-330 to be resolved by the ASM ( |Heinke et al 
|2001| l. Furthermore, M 15 X-2 is close to another source, AC 
211, in the globular cluster M 15. The ASM does not resolve 
these two sources, but because AC 211 is 2.5 times fainter than 
M 15 X-2 ( |Dieball et al.|2005| > we still use the combined X-ray 
data (known as 4U 2127+1 19) in this paper. 



2.3. Data analysis 

In order to calculate the luminosity distribution of a source, we 
use two methods, each based on individual dwells. The first 
method selects only those dwells that have a signal-to-noise ra- 
tio [S /AO above a fixed threshold, while the second method uses 
time-averaged luminosity to estimate the faint-end variability. 

A commonly used technique is grouping dwells into time 
bins of a fixed duration to statistically improve the S IN, thereby 
revealing fainter behavior. We do not apply this because in order 
to benefit from this, the bin duration should be on the order of 
days, and variability on shorter timescales will be lost. In partic- 
ular the bright-end behavior will be affected as high count rate 
dwells are combined with faint or noisy dwells. 



2.3.1. Signal-to-noise threshold 

The bright-end part of the luminosity distribution is best revealed 
by using the shortest available observations for maximum sensi- 



http : / /heasarc . gsf c . nasa . gov/ docs/ xte/asm_ 



products .html 



tivity, i.e., individual dwells. The downside of this is that many 
individual dwells have a count rate near the noise level. 

We use a strict criterion to select statistically significant 
dwells. For a given source, only dwells with S /N > cr are se- 
lected, where cr is found by assuming that a fraction l/n of the n 
dwells lies outside cr times the standard deviation of a Gaussian 
probability distribution (For bes et al.|2010 l 



i.e. 



1 - erf 



\ V2/ n 



(1) 



The factor 2 on the right-hand side appears due to asymmetry; 
only positive rates are considered. Dwells failing this criterion 
are discarded. The critical S /N ~ 4.2 for n = 7 X 10 4 , a typi- 
cal ASM value. The remaining sample is statistically expected 
to contain 1 data point that is not 'real' (in the sense that a dis- 
tribution fully dominated by noise would be expected to have 1 
outlier in its bright-end tail sufficiently far from the mean). This 
is the advantage of our method; the selected dataset consists of 
(almost) exclusively 'real' detections and therefore it can be used 
to find a lower limit of the system's average bolometric luminos- 
ity. However, this dataset is not complete; at the times when the 
detection has an insufficient S /N to pass the criterion, the actual 
luminosity of the source is unlikely to be zero. This luminosity 
can be estimated by an extrapolation of the known part of the 
distribution ( |Press et aXJ l992 chap. 3). We adopt the simple ap- 
proach of extrapolating the bright-end part towards lower lumi- 
nosities by a power-law function, as suggested by the relatively 
straight curves of the luminosity distribution when plotted as a 
log-log graph (see also Sect. 



3.1.1 



The interval over which we 
calculate the slope ranges from a fraction that is a factor 10 2 
below the fraction of time during which the source is detected 
significantly, to a fraction of 10~ 4 ' 5 . The faint-end and bright-end 
parts are excluded for practical reasons that will be explained in 
Sect ~ 



3.1 



From the extrapolation, we can integrate the (long- 
term) average luminosity. The validity of this extrapolation can 
be verified by comparing the resulting time-averaged luminosity 
with the time-averaged ASM luminosity and theoretical models. 

If the logarithmic slope of a cumulative luminosity distri- 
bution (i.e., a time-above-luminosity against luminosity curve 
is /?, then the luminosity distribution (i.e., time-at-luminosity 
against luminosity curve) function slope is (3 - 1 . Combined with 
the minimum luminosity L () , which is the luminosity L where the 
cumulative distribution has a fraction 1, and the Eddington lumi- 
nosity Lgdd (see Fig. [T|), this yields th e time-averaged bolometric 



luminosity (/?<-!) (Bremaud 1997, chap. 3) 



f LEM LxLP- l dL 

JLq 



JLr, 



T l3+\ _ T /3+l^ 
^Edd ^0 

^Edd ^0 ' 



jS+1 



La. (2) 



The approximation holds if {L^nl Lq)^ +1 <k 1 (which is the case 
for all known UCXBs except 2S 0918-549). 

2.3.2. Faint-end power-law function 

Instead of assuming that the source luminosity behavior can be 
described by a single power-law function as we did in Sect. 
|2.3.1[ here we assume that the source behavior at the times that 
the source is fainter than L slg (defined as the faint-end boundary 
of the range over which f3 is calculated) can be described by a 



2 Not to be confused with the X-ray luminosity function of a popula- 
tion of X-ray binaries. 
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Table 1. Known ultracompact X-ray binaries and candidates. The four groups separated by horizontal lines are confirmed UCXBs, 
UCXB candidates with tentative orbital period, UCXB candidates with a low ratio of optical to X-ray luminosity and UCXB 
candidates with persistent behavior at low luminosity, respectively. The columns list source name, orbital period P or b, distance d, 
location within the Galaxy, chemical composition and indication of accretor type ('msp' = millisecond pulsar, 'burst' = Type 1 
X-ray burst). Sources with coordinates very near the Galactic Bulge are assumed to be located in the Bulge. Distances are taken 
from Liu et al. ( 2QQ7\ unless indicated otherwise. In the references column, 'p' refers to orbital period, 'c' to composition, 'd' to 
distance and 'a' to accretor signal. 
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second power-law function, with slope y < 0, which in general 
is different from the power-law slope /3 between L s ; g and Ledd 
(see Fig. [Til. The two power-law functions meet in L s [ g . We solve 
for y andthe corresponding value of L ( * (the luminosity where 
the cumulative distribution based on both slopes /? and y has a 
fraction 1) by assuming that the integrated luminosity between 
L ( * and Ledd equals the time-averaged ASM bolometric luminos- 
ity Lasm (using all dwells), 



J^LxfiDdL 
£ M f(L)dL 



(3) 



where the luminosity distribution function 



f(L) = 



L y-\ 

Slg 



if L* < L < L sig , 
if L si „ <L< L Edd . 



(4) 



sl2 to find the factor Z7- 13 in front 

S1 S sig 



We used the continuity of / in L. 

A steep y-slope (very negative y) means that the source has 
a relatively constant luminosity when it is fainter than L s j g , and 
spends much time at luminosities not far below L s j g . A flat y- 
slope on the other hand points to a generally faint source that 



It follows that also the cumulative luminosity distribution is contin- 



uous in L„i 
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l/n 



Luminosity (logarithmic) 

Figure 1. Sketch of a cumulative luminosity distribution. The 
highly significant bright-end distribution (solid) and its linear ex- 
trapolation (dash-dotted) have slope f3 and are discussed in Sect. 



2.3.1 The extension of the bright-end distribution with slope y 
(dashed) is discussed in Sect. |2.3.2| The variables on the axes 
are defined in these two sections. In general the y-slope can be 
either flatter than (shown), equal to, or steeper than the /J-slope. 



spends a considerable amount of time far below L s i g , and is 
highly variable. 

This method does not work for two sources, XTE J 1807-294 
and XTE J0929-314, because L ASM < for thoseQ Hence, the 
contribution to Lasm from luminosities below would have 
to be negative. In reality this suggests that these sources must be 
very faint or off when not detected significantly, but more precise 
conclusions are not possible. 

2.4. Theoretical evolution 

Mass transfer in UCXBs is driven by angular momentum loss 
via gravitational wave radiation, since angular momentum loss 
via magnetic braking is much weaker at such a short orbital pe- 
riod (Landau & Lifs hitz|1975| Rappaport et al.|1983 i. Therefore 
the mass transfer rate strongly depends on the donor mass and 
orbital separation. During most of the evolution, the donor is de- 
generate and expands with mass loss. Since the donor keeps fill- 
ing its Roche lobe and the ratio between donor and accretor mass 
decreases, the orbital separation expands with time. The decreas- 
ing donor mass and increasing orbital separation cause the grav- 
itational wave emission and the corresponding mass transfer rate 
to decrease with time. 

A fit of the mass transfer rate M (defined > 0) in an 
UCXB with a zero-temperature helium white dwarf donor and 
an 1 .4 M Q neutron star accretor based on theoretical evolution- 



ary tracks (van Haaften et al. 2012b| l and optimized for an orbital 
period P OI \, < 1 hr is 



M = 1.2 x 10 



12 (Porb\ 

l~hr~i ' 



M Q yr 



(5) 



4 This is due to a systematic error. Faint sources, or sources far out- 
side the center of the field of view, have an increased probability to 
suffer from an imperfect solution for the brightness and position of the 
sources in the ASM field of view, and can end up with a consistently 
negative count rate, resulting in a negative average count rate over all 
dwells (Levine et al. 1996} . 



A simple estimate of the time-averaged bolometric luminosity L 
is 



L = 



GM,M 



2R a 



(6) 



where G is the gravitational constant and M a and R a the accretor 
mass and radius, respectively. For M a = 1 .4 M Q and R a = 12 km 
Eqs. |5]l and |6]l combine to 



L = 6 x 10 



33 



erg s 



(7) 



which is shown in Sect. 13.21 

Young UCXBs can contain a helium burning donor (Savonije 
et al. 1986[), in which case the mass transfer rate at a given orbital 



period is much higher than in the case of a white dwarf donor, 
due to the higher donor mass. During this evolutionary stage the 
orbit shrinks because the donor shrinks sufficiently rapidly as it 
loses mass. Once the donor has lost enough mass, core fusion is 
extinguished and a hot white dwarf donor remains. The system 
goes through a period minimum. Further mass loss and cool- 
ing cause the evolution to approach the zero-temperature white 
dwarf donor evolution. 

UCXBs may also form via an evolved main sequence donor 
scenario (Podsiadlowski et al. 20 Q2)>, but only fo r highly fine- 
tuned initial conditions ( |van der Sluys et al.|2005 1, therefore we 
do not consider this scenario. 



3. Results 

3. 1 . Signal-to-noise threshold 

By selecting only dwells that contain a highly significant source 
detection as described in Sect. 2.3. 1| a sensitive bright-end cu- 
mulative luminosity distribution can be constructed, as shown in 
Fig. [2] This figure shows how often a source emits above a given 
luminosity. Strictly speaking it shows the fraction of dwells, but 
this is similar to the fraction of time because of their unbiased 
timing. Here, the total number of dwells includes those that were 
discarded. 

From the vertical locations of the faint-end points of the 
curves it can be seen that typically only ~0.1 - 1% of the dwells 
pass the high-threshold criterion, emphasizing the fact that only 
the brightest phases of a source can be studied using this method. 
The short orbital period systems are generally visible more often 
than the long (in the context of ultracompact binaries) orbital pe- 
riod systems, as expected from theory. 4U 1820-30 (orbital pe- 
riod 1 1 .4 min) is the only source that is clearly observed nearly 
all the time (at 88% of the time). Most UCXBs are strongly vari- 
able, in the sense that curves are far from vertical. The most ex- 
treme case is 2S 0918-549 (17.4 min), which has a very flat 
logarithmic slope of /3 ~ - 1 in Fig. [2] Most slopes, however, are 
steeper with /3 < -2 (see Table|2]l. 

Since luminosity is used on the horizontal axis rather than 
flux or count rate, and given that the faintest data points of the 
cumulative distributions have a roughly similar count rate (deter- 
mined by the instrument's detection limit), the source distance 
determines at which luminosity the faint end of each curve is lo- 
cated. The UCXB candidate 4U 0614+09 1 is by far the neares t 
system at an estimated distance of 3.2 kpc (Kuulkers et al. 2010| l, 
therefore this source can be observed down to a lower intrinsic 
luminosity than the more distant sources. 

An artifact of the technique used in Fig. [2] can be seen at 
the faint ends of several curves, where at low luminosities the 



4 



van Haaften, Voss & Nelemans: Long-term luminosity behavior of 14 ultracompact X-ray binaries 



3.1.1. Faint-end behavior 
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Figure 2. Cumulative probability distribution of UCXB bolo- 
metric luminosity based on ASM data (fraction of time that a 
source emits above a given luminosity) for 14 UCXBs includ- 
ing two candidates with a tentative orbital period. Only count 
rates that are almost certainly real (their significant detection is 
almost certainly not the result of noise) have been used (Sects. 
|2.3.1| and |3.1[ ). The labels indicate the orbital period in min of 
the corresponding systems. Different line styles are used only to 
distinguish different systems. 



slope flattens towards a fraction of time usually much lower 
than 1 . This behavior most likely does not have a physical ori- 
gin and happens because at these faint luminosities, by chance 
a small number of data points will have an unusually low error 
and therefore pass the threshold, while the vast majority of data 
points with approximately the same count rate fails the thresh- 
old. This extends the curves further to low luminosities almost 
horizontally. This effect is more pronounced for sources that 
have more significant data points (corresponding to a fraction of 
time > 10~ 3 ), as those have a higher probability to have outliers 
(detections with an unusually low error). 

Even though the observation schedule of the ASM is prede- 
termined, irregular observation intervals may introduce an extra 
bias when an above- or below-average observation frequency co- 
incides with a period of abnormal source behavior. Since Fig. [2] 
uses individual dwells regardless of their timing, the shape and 
position of these distributions could potentially be different from 
their shape and position in the case of using bins of fixed dura- 
tion. This effect appears to be weak since there is almost no cor- 
relation between interval duration and count rate. At a fraction 
below ~ 1CT 4 the slopes of the curves become less reliable due 
to the low number of dwells the curves are based on, in particu- 
lar there is the possibility of atypical interval lengths surround- 
ing the highest count rates. Even in the case of regular sampling 
there would be some bias at the highest count rates, because even 
then small number statistics allow a source to behave differently 
during observations compared to in between observations. 

Although some luminosities in Fig. [2] appear to go above the 
Eddington limit for a neutron star, we note that given the uncer- 
tainties in both the distance and the ASM count to flux conver- 
sion, these luminosities are consistent with being sub-Eddington 
within errors. 



The nearly constant slopes over a large part of the curves in 
Fig. [2] (excluding the faint-end flattened part and the statisti- 
cally uncertain bright-end part of these curves for the reasons 
discussed in Sect. 



3.1 



suggest the possibility that this behavior 
may continue into the faint-end phases of the sources that have 
not been detected reliably. By extrapolating the middle parts of 
these curves with a power law, the time-averaged luminosity can 
be expressed in essentially only the slope and the minimum lu- 
minosity, unless the slope is close to /3 — -1 (Sect. 2.3.1 1. In 
Sect. 3.2 the result will be compared with the time-averaged lu- 
minosity and with theory. 

3.2. Average luminosity 



The variability behavior of the UCXBs presented in Sect. 3.1 can 
be summarized and compared by using the time-averaged lumi- 
nosity. Figure [3] shows several bolometric luminosity estimates 
for the known UCXBs. The average luminosity over all dwells 
(circles) can be seen as the best estimate of these three taking all 
uncertainties into account, unless, of course, when the average 
count rate is negative. In general, the luminosity estimates based 
on power-law extrapolation (Sect. 3 - 1 - 1 1 > and the time-averaged 
value over all dwells are of the same order of magnitude. 

Figure [3] also compares the average-luminosity estimates 
with theoretical UCXB tracks, with both a zero-temperature 
white dwarf donor and a helium burning donor. The helium burn- 
ing donor system has typical initial conditions among the ranges 
simulated by Yungelson (2008]). 

The average luminosity of the short-period UCXBs (below 
~ 30 min) agrees well with the theoretical model based on the 
mass transfer rate (Of these, M 15 X-2 is contaminated by a 
fainter source, see Sect. |2.2.1[ ). Based on the ASM data, there 
is no reason to suspect variability on timescales longer than 16 
yr in these sources. The extrapolated power-law estimate closely 
matches both the average luminosity and the theoretical predic- 
tion. 4U 1820-30 (1 1.4 min) is the only source with a very well- 
determined average luminosity, which fits well with the repre- 
sentative helium burning donor track. As this source may have 
a negative period derivative ( Tan et al. |1991 1, it could corre- 
spond to the upper branch of a helium donor track, but this is 
uncertain. 2S 0918-549 (17.4 min), which has the flattest slope 
in Fig. [2] has a correspondingly low luminosity estimate based 
on the extrapolated power law, much lower than even the zero- 
temperature white dwarf donor model. This means that either 
the source has gone through an unusually faint period of time 
the last 16 years, or the extrapolation underestimates its faint- 
end luminosity. 

The vertical dotted line in Fig. [3] represents the critical or- 
bital period for thermal-viscous disk instability, in the case 
of a |Dunkel et al.| ( |2006| l helium accretion disk and a zero- 
temperature white dwarf donor. It lies in a large apparent period 
gap ranging from 23 to 40 min. 

Above this gap, the four UCXBs with orbital periods be- 
tween 40 and 45 min have lower bounds (triangles) close to the 
zero-temperature donor track. This is consistent with them be- 
ing (very) faint when not detected, i.e., transient, as is expected 
for systems with these orbital periods. The power-law extrapo- 
lations (stars), which give a much higher luminosity estimate, 
would not apply then. 

However, the three UCXBs, including one candidate, with 
an orbital period longer than ~ 50 min all have a lower bound 
on their luminosity that exceeds both the zero-temperature white 
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Figure 3. UCXB time-averaged bolometric luminosity against 
orbital period. The solid curve represents the modeled evolution 
of an UCXB with a zero-temperature helium white dwarf donor 
(Eq. [71. The dashed curve represents the modeled evolution of 
an initially detached system with an 0.65 M Q helium burning 
donor starting at an orbital period of 100 min, which starts mass 
transfer via Roche-lobe overflow at an orbital period of 46 min 
( Yungelson 2008|l. Both systems have an 1.4 M Q neutron star ac- 



cretor. The triangles represent the absolute lower bounds on the 
luminosity based on only the highly significant ASM dwells and 



assuming zero luminosity at all other times (Sects. 2.3.1 and 3.1 
Fig. [2j. The stars above these represent the power-law extrapo 
lation of the bri ght-end part of the cumulative luminosity distri- 
bution La (Sect. 



3.1.1 



The circles represent the time-averaged 
bolometric luminosities based on all dwells Zasm- Filled sym- 
bols correspond to confirmed UCXBs, open symbols correspond 
to UCXB candidates with a tentative orbital period. Symbols at 
the same orbital period correspond to the same source and are 
connected by a gray line for clarity. For two sources the aver- 
age count rate is negative so the circle is not shown. The ver- 
tical dotted line at 28 min represents the orbital period above 
which a helium accretion disk becomes subject to the thermal- 
viscous instability, based on the zero-temperature white dwarf 
donor model by van Haaften et al. ( 2012b| l. 



dwarf track and the (extrapolated) degenerate part of the helium 
donor track. Different initial conditions can result in a slightly 
higher helium donor track (Yungelson 2008) so some of these 
systems can still be consistent with having a (partially) degen- 
erate donor. One source, 4U 0614+091, cannot be explained 
by having a semi-degenerate donor as its lower bound stands 
out above the degenerate track too much. 4U 0614+091 has an 
uncertain orbital period, but if this period is correct then this 
suggests the system has a helium burning donor. Another pos- 
sibility is that this source was exceptionally bright during the 
past 16 yr, and that its average luminosity over a longer time 
is lower. Similarly, the time-averaged ASM luminosities of all 
three UCXBs with an orbital period above ~ 50 min, as well as 
that of XTE J1751-305 at 42 min indicate that these systems are 
much brighter than predicted by the theory for UCXBs with a de- 
generate donor, but still fainter than expected for UCXBs with a 
helium burning donor. In several cases the time-averaged lumi- 
nosity is nearly equal to the power-law extrapolation estimate, 
similar to the short-period systems. It is unlikely that all three 
known systems with orbital period > 50 min are in a bright state 



Table 2. Power-law slopes /3 and y describing the cumulative lu- 
minosity distribution (see Sects. 2.3.1 and 2.3.2 1 for the known 
UCXB population based on the RXTE ASM data. The two sys- 
tems below the line are UCXB candidates, with an uncertain or- 
bital period. If no y-value is given, the source is too faint on 
average for this method to work. 



Ntimc 


P u (mm) 

1 orb V 111111 / 


P 


y 


4U 1820-30 


11.42 


-1A 


-3.7 


4U 0513-40 


17 


-2.0 


-2.3 


2S 0918-549 


17.4 


-1.03 


-2.5 


4U 1543-624 


18.2 


-3.8 


-4.3 


4U 1850-087 


20.6 


-1.9 


-2.1 


M 15 X-2 


22.58 


-1.6 


-3.0 


XTE J 1807-294 


40.07 


-2.4 


N/A 


4U 1626-67 


41.4 


-3.0 


-2.5 


XTE J175 1-305 


42.42 


-2.5 


-0.55 


XTE J0929-314 


43.58 


-2.8 


N/A 


4U 1916-05 


49.48 


-2.8 


-2.5 


SWIFT J 1756.9-2508 


54.70 


-1.3 


-1.7 


4U 1728-34 


10.77? 


-6.3 


-0.66 


4U 0614+091 


51.3? 


-3.2 


-4.1 



by chance. An explanation could be that the theory described in 



Sect. 2.4 underestimates the time-averaged mass transfer rate of 
long-period systems. Alternatively, UCXBs at long periods may 
be variable on a timescale much longer than 16 yr, so that during 
the ASM observations some systems were much brighter and 
many others (much) fainter than the average luminosity. Such 
a selection effect would lead us to observe only the atypically 
bright systems. 

3.3. Faint-end power-law function 



2.3.1 stars in 



If the extrapolated power-law method Lp (Sect. 
Fig. [3} underestimates the time-averaged luminosity LAsiv|^](cir- 
cles in Pig. [3j, then the slope of the cumulative luminosity distri- 
bution is expected to steepen below L s j g , i.e., y < ft. This means 
that the source is brighter and more constant at low luminosi- 
ties (fainter than the significant detections) than suggested by 
extrapolating the behavior during the bright phases. Conversely, 
if the extrapolated power-law method overestimates the time- 
averaged luminosity then y > ft and the slope flattens below L s i g . 
In this case the source is more variable when faint than when 
bright. Table [2] shows the values of /3 for the known UCXBs 
as well as y for the sources with a sufficiently high (in prac- 
tice, positive) time-averaged luminosity. Steepening of the faint- 
end power-law slope happens for all short-period sources (be- 
low ~ 30 min) except 4U 1820-30 at 11.4 min. The latter may 
be caused by the fact that this source is detected significantly a 
very high fraction (88%) of the time. On the other hand, almost 
all long-period sources (above ~ 30 min) behave the opposite 
way, with a flattened faint-end slope. The one certain exception 
is SWIFT J 1756.9-2508 at 54.7 min, while 4U 0614+091 would 
be a second exception if its orbital period would be confirmed. 



4. Implications for population studies 

The results concerning the luminosity behavior presented in 
Sect.[3]can be used in modeling a population of UCXBs. Based 



5 The time-averaged luminosity is calculated based on the average 
count rate over all individual dwells (regardless of their timing), which 
is a good approximation given the unbiased timing and large datasets. 
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on its evolutionary state, a modeled UCXB can be associated 
with a real UCXB, and its luminosity and variability can be 
estimated empirically by drawing from the observational lumi- 
nosity distribution, without need to understand the theory be- 
hind the variability. This method rests on the assumption that the 
available observational record is typical for such systems, which 
may not be the case for the long-period UCXBs observed by the 
ASM, as shown in Fig.[3]and noted in Sect. |3.2| 

Since many UCXBs are only visible a fraction of the time 
above a given luminosity as shown by Fig. [2] at any given time 
we will only see a fraction of the UCXB population, using an in- 
strument with a given sensitivity. If the instrument is the RXTE 
ASM itself, only a small fraction of the population can be seen 
at any given time. For a more sensitive instrument, the shape 
of the luminosity distribution determines the characteristics of 
the visible population. A flat slope at the faint end of the cu- 
mulative luminosity distribution implies a large, faint popula- 
tion, while a steep faint end implies a small, relatively bright 
population. Theoretical models relatively accurately predict the 
integrated amount of energy that is emitted by a particular pop- 
ulation (based on time-averaged mass transfer rate and gravita- 
tional energy release), but are less accurate when it concerns the 
variability of sources, which relates mainly to complex accretion 
disk behavior. 



5. Discussion and conclusions 

We have investigated the long-term X-ray light curves collected 
by the RXTE ASM for 14 ultracompact X-ray binaries including 
two candidates. 

All UCXBs have a significantly varying flux during their 
brightest phases (above 10 37 erg s _1 ) as shown by the rela- 
tively flat slopes in Fig. [2] By comparing their (integrated) time- 
averaged bolometric luminosities with theoretical estimates (Fig. 
[3]), we conclude that for short-period systems (systems with 
an orbital period below ~ 30 min) these flat slopes proba- 
bly continue at more typical luminosities (below 10 37 erg s _1 ), 
since these relatively flat cumulative luminosity distributions are 
needed to match the observed and theoretical average luminosi- 
ties. Longer-period systems must have a much flatter slope than 
short-period systems at luminosities below 10 37 erg s _1 in order 
to match the theoretical models, in other words, they are more 
variable when faint relative to short-period systems. This follows 
from the steep decrease in theoretical luminosity with increasing 
orbital period. In general, all sources probably vary strongly at 
all luminosities. This confirms that long-term observations are 
indeed needed in order to compare a source's luminosity to a 
theoretical model of the mass transfer rate. 

Considering the time-averaged ASM luminosity over all 
dwells, for the long-period systems a single power law yields 
an average luminosity that is too high. Instead, the slope must 
become flatter at lower luminosities, as already suggested by the 
theoretical model. The short-period systems on the other hand 
require a steeper low-luminosity slope. 

Furthermore, several long-period systems have an average 
luminosity that is much higher than predicted by theory in the 
case of degenerate donors. A possible (partial) explanation is 
that the low-mass donors in these UCXBs are being heated and 
evaporated by irradiation from the millisecond pulsar they are 
orbiting. The angular momentum loss from the system resulting 
from a stellar wind from the donor enhances the mass transfer 
and hence increases the luminosity, as has been demonstrated 
in explaining the millisecond pulsar binary system PSR J 1719- 



1438 ( |van Haaften et al.|2012a[ ). A larger donor radius in itself 
also leads to a higher mass transfer rate. 

There is no clear distinction between transient and persis- 
tent behavior in the sense that the estimated average luminosities 
are similar for sources over a significant range in orbital peri- 
ods, which is partially a result of the similar bright-end slopes in 
the cumulative luminosity distributions. However, the trend that 
long-period systems typically have a flatter slope when they are 
faint compared to when they are bright can be seen as transient 
behavior. The variability of short-period UCXBs is more consis- 
tent with a single power-law cumulative luminosity distribution 
that holds up to a fraction of time of 1, which can be understood 
as 'persistent' behavior. 

Some UCXBs seem too bright to have a degenerate donor. 
Especially for the long-period systems, orbital period derivatives 
would be very useful to decide whether the donor is helium burn- 
ing or degenerate. This knowledge in turn would limit the range 
of theoretically predicted time-averaged luminosities, and thus 
constrain the faint-end luminosity distribution. 

The orbital period of 4U 0614+091 may be much shorter 
than the ~ 50 min that is suggested by several observations (see 
Table[T i. Both the high lower bound on the luminosity compared 
to the degenerate-donor model for long orbital period UCXBs, 
and the apparent steepening of the slope of the cumulative lu- 
minosity distribution at low luminosities, resemble the behavior 
of short-period systems. At least the latter argument is not con- 
clusive, given that SWIFT J1756. 9-2508 combines a 54.7 min 
orbital period with a steepening slope at low luminosities. 

Even though the RXTE ASM observations are the longest 
available, there may very well be variability on timescales 
(much) longer than 16 yr, which could imply that the data for 
some of these systems available today are atypical. This could 
affect the average luminosity as derived by each of the meth- 
ods, including the lower bounds. As for the long-period UCXBs, 
apart from being consistently brighter than theoretical estimates 
for systems with degenerate donors, variability on very long 
timescales (hundreds or thousands of years) would ensure that 
during any 16 yr period, a small part of the population is excep- 
tionally bright. This could be the population we observe. 

The Monitor of All-sky X-ray Image Gas Slit Camera (MAXI 
GSC) has a higher stated sensitivity (Matsuoka et al. 2009) than 
the ASM over a shorter period of time of 2.5 yr. Because of the 
much longer observation baseline, we decided to use only the 
ASM data. We compared the MAXI data with the ASM data and 
found that the MAXI sensitivity is better than the ASM data for 
7 UCXBs with known (or tentative) orbital period, similar for 2, 
worse for another 2 and unavailable for 4 UCXBs. Overall the 
difference was not large enough to outweigh the shorter obser- 
vation baseline. 

We will use the luminosity distributions found in this paper 
in a forthcoming study of the UCXB population in the Galactic 
Bulge, to estimate the X-ray luminosities of systems predicted 
by a population model. In the context of predicting the observ- 
able population of UCXBs, the Galactic Bulge Survey ( |Jonker| 
et al.pOlT) may constrain the number of faint sources, which 
also gives information on the luminosity distribution. 
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